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Q ■ ABSTRACT 

\ We combine in a self-consistent way the constraints from both gravitational lensing and stel- 

i lar kinematics to perform a detailed investigation of the internal mass distribution, amount of 

' dark matter, and dynamical structure of the sixteen early-type lens galaxies from the SLACS 

Survey, at z = 0.08 - 0.33, for which both HST/ ACS and NICMOS high-resolution imaging 
' and VLT VIMOS integral-field spectroscopy are available. Based on this data set, we analyze 

' the inner regions of the galaxies, i.e. typically within one (three-dimensional) effective radius 

re, under the assumption of axial symmetry and by constructing dynamical models supported 
\^ • by two-integral stellar distribution functions (DFs). For all systems, the total mass density dis- 

fSJ I tribution is found to be well approximated by a simple power-law ptot '« (with m being 

CNJ . the ellipsoidal radius): this profile is on average slightly super-isothermal, with a logarithmic 

■ slope (/) - 2.074+QQ4J (errors indicate the 68% CL) and an intrinsic scatter cry - 0.144+qqj4, 

[ and is fairly round, with an average axial ratio (q) = 0.77 + 0.04. The lower limit for the dark 

matter fraction (/dm) inside rg ranges, in individual systems, from nearly zero to almost a 
half, with a median value of 12%. By including stellar masses derived from stellar popula- 
tion synthesis (SPS) models with a Salpeter initial mass function (IMF), we obtain an average 
. ^ /dm = 31%, and the corresponding stellar profiles are physically acceptable, with the excep- 

I tion of two cases where they only marginally exceed the total mass profile. The /dm rises to 

^ , 6 1 % if, instead, a Chabrier IMF is assumed. For both IMFs, the dark matter fraction increases 

_ with the total mass of the galaxy (correlation significant at the 3-sigma level). Based on the 

intrinsic angular momentum parameter calculated from our models, we find that the galaxies 
can be divided into two dynamically distinct groups, which are shown to correspond to the 
usual classes of the (observationally defined) slow and fast rotators. Overall, the SLACS sys- 
tems are structurally and dynamically very similar to their nearby counterparts, indicating that 
the inner regions of early-type galaxies have undergone little, if any, evolution since redshift 
z«0.35. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: kinematics and dynamics — 
galaxies: structure — gravitational lensing 



1 INTRODUCTION 

Early-type galaxies are among the fundamental constituents of the 
local observable Universe, contributing to most of the t otal stellar 
mass jpukugita. Hogan. & PeebTe^l 19981 : lRenzini|[200^) . and it is 
therefore only natural that these systems have become the focus of 

* E-mail: mbamabe@stanford.edu 



a vast amount of observational and theoretical studies during recent 
years. 

Understanding the formation and evolution mechanisms that 
can generate such a remarkably homogeneous class of objects re- 
mains among the most important open questions in present-day as- 
trophysics: although the standard scenario according to which mas- 
sive elli pticals are built-up via hierarchical merging of disk gala xies 
(see e.g |Toomrd[T9T7llWhite&Frenklll99ll : ICole etal.ll2000l) has 
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been steadily gaining consensus throughout the last decades, only 
recently the progress in high-resolution numerical simulations has 
allowed to start investigating i n detail th e internal structural prop- 
erties of realistic systems (e.g. 



1003, Naabetal. 2007, 



jMeza et al . 

Jesseit et al.l2007l.lOfiorbe et al.l20()7llGon zalez-Gar cia et al.l200S 
Johansson. Naab. & Ostrikeill2009llLackner & OstrikeilbOlol) 



Clearly, in view of these developments, the availability of 
stringent observational constraints has become even more cru- 
cial. Luckily, at least as far as nearby ellipticals are concerned, 
the observational and modelling effort has been very intense 
in the past few years. Stellar dyna mics represents the most 
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ring (e.g. lOosterloo et al. 2002; Weijmans et al. 2008), and hot X- 
ray gas emission (e.g^Matsushita et al. 1998; Fukazawa et al . 20061 ; 
iHumphrev & Buotell20rol ; lOas et alj|2010l ; IChurazov et al.ll2010l) . 
In general, from these analyses there is mounting evidence that 
early-type galaxies are embedded in dark haloes, whose contri- 
bution is often significant already inside one effective radius, and 
that their total density profiles are roughly isother mal within the in- 
ner regions (but not nece ssarily at larger radii, see lde Lorenzi et afl 
l2009llDutton et al.ll2010l) . 

Ideally, if such studies could be extended to elliptical galax- 
ies beyond the local Universe, it would become possible to pro- 
vide strong constraints also to the intermediate steps of numeri- 
cal simulations throughout redshift, rather than just to their z = 
end-products. The current observational limitations, unfortu- 
nately, make it prohibitively difficult to apply the traditional tech- 
niques to distant systems, i.e. at z x 0.1 and beyond. In par- 
ticular, the impossibility to reliably measure the higher-order ve- 
locity mo ments needed to disentangle the mass-ani sotropy de- 
generacy ( lGerhard|[l993l;lvan der Marel & Franxll 19931) constitutes 
a serious hindrance for dynamical studies. On the other hand, 
however, at higher redshifts strong gravitati onal lensing (see e.g. 
ISchneider. Kochanek, & Wambsgansi l200d) may become avail- 
able as an additional probe, opening up new possibilities to investi- 
gate the mass distribution of galaxies. 

The main obstacle with this approach, namely the diffi- 
culty in discovering a suitable number of early-type galaxies 
that also act as strong gravitational lenses, has been rem edied 
by the Sloan Len s ACS Survey, SL ACS teolton et alj 
2008a, 2008b, Koopmans et alj I 



Gavazzi et all 120071 , |2008| , lAuger 



200^. Ft 
er et al.l I 



T reu et al. 
20091 




which 



has yielded a large sample of almost a hundred such objects. 
As highlighted by these studies, as well a s by a number of 
other works (see Treu & Ko opmans' 2002', 'Koopman s & Treiil 
2OOI iTreu & Koopmani [2004, R usin & Kochane j 1 20051, 
Jiang & Kochanekl I2OOI iGrillo et alj |2008. Koop mans et alj 
20091 iTrott et alj|2010l Ivan de Ven et ^l20ia iDutton et alj|20rir 



Spiniello et alj 2011 ), the combination of gravitational lensing 
and stellar dynamics provides a particularly powerful and robust 



method to determine the total mass density distribution of distant 
galaxies. 

A potential limitation in the joint analysis of the SLACS 
galaxies mentioned above is that all the information about the dy- 
namics is extracted from a single data point, i.e. the average stel- 
lar velocity dispersion measured within a 3 arcsec diameter aper- 
ture, obtained from SDSS spectra. Moreover, the method is not 
fully self-consistent, in the sense that axial symmetry is assumed 
for the lensing modelling, but not for the dynamical modelling, 
which is based on spherical Jeans equations. In order to address 
the concerns that these approximations might lead to biased re- 
sults, the SLACS data set has been complemented, for about 30 
lenses, with extended, two-dimensional kinematic information, ob- 
tained either with VLT VIMOS integral field unit (IFU) or from 
LRIS Keck long-slit spectroscopic observations (with the slit posi- 
tions offset along the galaxy minor axis in order to mimic integral 
field capabilities). In addition, we have developed a more rigorous, 
self-consistent modelling technique (the cauldron code, based on 
the superposition of two-integral distribution functions, see Sect.|3] 
and lBamabe & Koopmansir2007l '), aimed at taking full advantage of 
the available data sets for each galaxy, namely the gravitationally 
lensed image, the surface brightness distribution and the velocity 
moments maps. This makes it possible to recover more accurate 
constraints on the slope and the flattening of the total mass den- 
sity distribution, to derive information on the dark matter fraction 
within the inner regions and to obtain insights on the global and 
local dynamical status of the system, so that this analysis consti- 
tutes, to a good extent, the higher redshift analogue of the detailed 
studies that are routinely carried out on local ellipticals. We have 
applied this method to half a dozen SLACS galaxies w ith VIMOS 
observations jCzoske et al.| l2008; Barnabe et al.' 2009ah , as well as 
to one of the Keck systems (Barnabe et al. 201Q). In this paper, we 
extend this in-depth analysis to the full sample of SLACS early- 
type lenses with available VIMOS IFU spectroscopy (we refer to 
Czoske et al. 201 1, in preparation, for a comprehensive description 
of the complete data set), a total of 16 systems at redshift z = 0.08 
to 0.35, and covering a wide range in mass and angular momen- 
tum. We also make use of stellar masses derived from SPS models 
jAuger et alj|2009l) to impose further constraints on the luminous 
distribution and estimate the contribution of the dark matter com- 
ponent. 

This paper is organized as follows; in Section |2] we provide 
an overview of the available data-sets. In Section |3] after recall- 
ing the basic features of the method for the combined analysis of 
the adopted model, we present the main results of the study of the 
SLACS lenses mass density profiles. The mass structure of the an- 
alyzed galaxies, in terms of luminous and dark matter components, 
is described in Section|4] while Section |5]deals with the recovered 
dynamical structure of the systems. We discuss our findings in Sec- 
tion|6]and we briefly summarize the main conclusions in Section|7] 
Throughout this paper we adopt a concordance ACDM model de- 
scribed by Dm = 0.3, £1^ = 0.7 and Hq = 100 A km s"' Mpc"' with 
h = 0.7, unless stated otherwise. 



2 OBSERVATIONS 

The joint lensing/dynamics analysis requires three types of input 
data; (i) high-resolution images are used to trace the surface bright- 
ness of the gravitationally lensed background galaxies; (ii) near- 
infrared images provide the surface brightness distribution of the 
stars in the lens galaxies; (iii) integral-field spectroscopic observa- 
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tions are used to derive two-dimensional maps of systematic veloc- 
ity and velocity dispersion of the stars in the lens galaxies. 

2.1 Imaging 

High-resolution images of the lens systems in the sample were ob- 
tained with Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). We use images taken through the F814W 
filter. Deep images are available for eight of the systems, while 
for the remaining nine only single-exposure snapshot images are 
available. The pixel scale of the images is 0.05 arcsec. The struc- 
ture of the gravitationally lensed background images was isolated 
by subtracti ng elliptical B-spline models of the lens galaxies' light 
distribution jSolton et aLll2008iJ) . 

NICMOS F160W images were used to obtain the surface 
brightness of the lens galaxies. In order to serve as luminosity 
weights to the kinematic maps, the images were matched in reso- 
lution to the spectroscopic observations (0.8 arcsec) and resampled 
from the original pixel scale of 0.05 arcsec to the pixel grid of the 
kinematic maps (0.67 arcsec). For three systems no NICMOS ob- 
servations were available. In these cases, the lens surface brightness 
was obtained from the B-spIine fits to the F814W images. 

2.2 Spectroscopy 

Two-dimensional spectroscopy of seventeen systems was obtained 
with the integral-field unit of VIMOS on the VLlQ Three sys- 
tems (J0037, J0912 and J2321) were observed with the HR-Blue 
grism(4000 - 6200 Aj^ , while the remainder were observed with 
the HR-Orange grism (5000 - 7500 Afl The field of view of 
27" X 27" was sampled with 40 x 40 spatial elements with a scale 
of 0.67 arcsec, the spectral resolution is R x 2500 (coixesponding 
to a velocity resolution of Ad k 90 - llOkms"' FWHM in the 
rest frame of the lens galaxies). Two-dimensional maps of system- 
atic velocity (with respect to the mean redshift of the lens galaxy) 
and velocity dispersion were measured by fitting stellar templates 
convolved with Gaussian line-of-sight velocity distributions to the 
spectra from spatial elements with a sufficient signal-to-noise ratio 
(S/N > 8 per spectral element of 0.644 A). 

A full description of the sample, observations, data reduction 
and the resulting kinematic maps is given in a companion paper 
(Czoske et al. 201 1, in preparation). 



3 COMBINED LENSING AND DYNAMICS ANALYSIS 

In this Section — after recalling the adopted modelling assumptions 
and the salient features of the method employed to self-consistently 
combine lensing and dynamics constraints — we present the core 
results of the study of the full sample of sixteen SLACS early-type 
lens galaxies for which data sets including both high-resolution 
imaging and VIMOS integral field spectroscopy are available. 

The theoretical framework of the joint lensing and dynamics 
analysis and the implementation of the method (i.e., the cauldron 

' One of the lens systems observed with VIMOS IFU, SDSS J1250-0135 
(identified as a late-type galaxy with morphological type SO/SA), has been 
excluded from the sample examined in this work due to its spiral arms par- 
tially appearing in the lens-subtracted image, severely hindering an accurate 
lensing analysis. 

2 ESO programme 075.B-0226, P.I.: L. Koopmans 
' ESO programme i77.B-0682, P.I.: L. Koopmans 



code) are described in full detail in ISamabe & KoopmansI ( l2007l 
hereafter BK07), to which we refer the interested reader. 



3.1 Method overview 

We describe the total mass density distribution of the lens galaxy as 
p(x; J]), where the vector x denotes the spatial coordinates and j] is 
a set of parameters that characterizes the density profile. The cor- 
responding total gravitational potential, (^(x), is calculated via the 
Poisson equation and utilized for both the gravitational lensing and 
the stellar dynamics modelling of the data set . By employing a pix- 
elated source reconstruction meth od (see e.g. IWarren & Dvell2003l . 
lKooDmansll2005l . lSuvu et alj2006l) for the lensing a nd some flavour 
of the Schwarzschild orbit superposition method dSchwarzschildl 
1 19791 : seeljhomas 20 1^ for a review of the developments of this 
technique) for the dynamics, both these modelling problems can be 
formally written in an analogous way as a set of regularized linear 
equatio ns, for which robu st solving techniques are readily available 
(see e.g. lPress et al.|[l992h . 

Clearly, for any given set of observations, the ultimate goal 
of the analysis is to determine the set of parameters that generate 
the 'best' density model, i.e. the most plausible in an Occam's razor 
sense. This can be achieved by embedding the linear optimization 
scheme within the framework of Bayesian inference, which pro- 
vides an objective way to quantify the probab ility of each model 
given the data (see e.g. lMacKa\il I999I . boOSl) . The best model is 
thus the maximum a posteriori (MAP) model, i.e. the one that max- 
imizes the joint posterior probability density function (PDF). 

While this method is very general, allowing in principle for 
any density distribution, in its current implementation — the caul- 
dron algorithm — we make further assumptions in order to achieve 
the computational efficiency that is needed for a thorough explo- 
ration of the parameter space. Within the cauldron code, there- 
fore, galaxies are modelled as axially symmetric systems, i.e. 
p(R, z), whose dynamics are described by a two-integral distribu- 
tion function (DF) / = /(£, L-) which depends on the two clas- 
sical integrals of motions, the energy E and the angular momen- 
tum along the rotation axis L-. This makes it possible to con- 
struct dynamical models in a matter of seconds, by means of a 
fast Monte Carlo numerical implementation developed by BK07 
of the two-integral S chw arzschild method originally in troduced by 
ICretton et al. I j 19991) and lVerolme & de Zeeuwl (12002!), which em- 
ploys two-integral components (TICs), rather than the commonly 
used stellar orbits, as dynamical building blocks. TICs can be 
viewed as elementary stellar systems completely specified by a par- 
ticular choice of E and L-; they are characterized by analytic radial 
density distributions and (unprojected) velocity moments, which 
makes them convenient and (compared to orbits) computationally 
inexpensive building blocks. 

Remarkably, the applicability of the code is not too severely 
limited by the assu mptions mentioned above: as shown in 
iBamabe et al. I l l2009ah . cauldron has been successfully tested for 
robustness by applying it to a complex system obtained from a nu- 
merical N-body simulation of a merger process, which obviously 
departs from the condition of axisymmetry. Nevertheless, several 
important global properties of the system (particularly the total den- 
sity slope and, when rotation is visible in the kinematics maps, the 
angular momentum and the ratio of ordered to random motions) 
are recovered quite reliably, typically within 10 to 25 per cent of 
the true value. 
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3.2 The galaxy model 



We model the total mass density profile of the lens galaxy with an 
axially symmetric power-law distribution, 

P(m) = (1) 

where po is a density scale, y' is the logarithmic slope of the profile, 
and the ellipsoidal radius m is defined as 



if2 



(2) 



here ag denotes the arbitrary length-scale and q is the axial ratio 
(q = I for a spherical galaxy, while oblate systems have ^ q < I). 
The c orresponding gr a vitatio nal potential is obtained by apply- 
ing the IChandrasekhaJ ( 1 19691) formula for homoeoidal density dis- 
tributions and can be written as a rather simple expression (cf. 
[Barnabe et al. 2009a) that just requires the computation of a one- 
dimensional integral. 

Despite its simplicity, the power-law model (typically with a 
slope close to the isothermal case, i.e. y' = 2) provides a good 
description of the total mass distribution of early-type galaxies 
over a large radial r ange, as reported by previous analyses of 



the SLACS system s i Koopmans et al 2 006|: Gavazzi et al 



Czoske et all (20081: iBamabe et al.ll2009d : iKoopmans et al 



2007 



200S 



Auger etal.ll2qi0ah. as well as a number of stell ar dynamics (e.i ^ 
Gerha rd et al. '200l|), gravitational len sing (e.g. Dye et al. 20081; 



Tu et al. 2009) and X-ray studies (e.g. iHumphrev & Bu ote 20 IQ) 



of ellipticals. 

Three physical parameters (the non-linear parameters i], in the 
notation of Sect. 13. U characterize the power-law model: the loga- 
rithmic density slope y', the axial ratio q and the lens strength Oq, 
a dimensionless quantity directly related to the normalization of 
the potential well (see Appendix B of BK07). Furthermore, four 
'geometrical' parameters are needed to describe the observed con- 
figuration of the lens galaxy in the sky, i.e. the inclination angle ; 
(with i = 90° being an edge-on system), the position angle and the 
sky-coordinates of the centre of the system. Luckily, the last three 
parameters are well constrained by the surface brightness distribu- 
tion of the lens galaxy and of the lensed images. Therefore, after 
determining their values by means of fast preliminary explorations, 
we maintain them fixed during the computationally expensive opti- 
mization and error analysis runs. 

The level of regularization, which controls the smoothness in 
the reconstructed background source and TIC weight maps, is set 
through three additio nal non-linear par ameters, usually known as 
hyper-parameters (see lSuvu et al.ll2006l and BK07 for a more tech- 
nical explanation). The optimal values of the hyper-parameters are 
also found via maximization of the joint posterior PDF. 



3.3 Bayesian inference and uncertainties 

For our analysis we follow the standard Bayesian framework de- 
scribed by, e.g., MacKay (2003) involving multiple levels of infer- 
ence: here we summarize the main steps involved in this approach, 
referring the interested reader to BK07 for a detailed explanation 
of the method. 

Let us denote the combined data set as d, the noise in the 
data as n and the considered hypothesis (i.e. the adopted model) 
as "Hixi), where the non-linear parameters // are the physical and 
geometrical quantities (such as the density slope, axial ratio, lens 
strength and inclination angle) that characterize the model. At the 
first level of inference, the model is assumed to be true (i.e. we fix 



one choice for the set of q parameters) and we aim to solve for w 
the linear problem 



d = Aw + n , 



(3) 



where w are the linear parameters — namely, the surface brightness 
distribution of the lensed background source and the weights of the 
dynamical building blocks — being mapped onto the observables 
by the operator A = A['7Y(7/)]. The construction of the mapping 
operator A for both the lensing and dynamics modelling constitutes 
the core of the cauldron code, and is described in BK07. From 
Bayes' theorem, the posterior probability of the parameters w given 
the data is written as 

Vv(d\w,q) Pr(a)|;u) 



(4) 



where the probability Pr(rf | w, if) is the likelihood and Pr(rf | w, ifi 
represents the prior. Here, we adopt as prior a curvature regular- 
ization function that describes the degree of smoothness that we 
expect to find in the solution; the amount of regularization to be ap- 
plied is set by the value of the so-called 'hyper-parameter' ji. The 
evidence Zi = Pr(rf|yu,//) = ^Y't(d\w,q)Vr{w\y)Aw is simply a 
normalization constant in Eq. (|4j, but it becomes important at the 
upper level of inference, where it appears as the likelihood. The set 
of parameters i«map that maximizes the posterior can be obtained 
by means of standard linear optimization techniques. 

Analogously, at the second level of inference the parameters 
are still maintained fixed and the posterior probability of the hyper- 
parameter is 

Pr(rf|;u,;7)Pr0u) 



Vv{n\d,Tj)- 



(5) 



here we adopt a scale invariant prior Pr(^) oc 1/yu (uniform in log) 
to formalize our ignorance of the order of magnitude of the hyper- 
parameter. 

Finally, at the third level of inference we can compare the dif- 
ferent models by studying the posterior probability distribution of 
the set of non-linear parameters 77, 

Pr(rf|/,)Pr(/7) 



Pr(77|rf): 



(6) 



and infer the parameters j/map for which the posterior is maxi- 
mized. For each galaxy, the corresponding MAP model consti- 
tutes the 'best model', meaning the most plausible joint set of 
model parameters given the data, and therefore throughout the pa- 
per we present the reconstructed observables, background source, 
TIC weights maps and inferred quantities that are obtained with the 
choice ri = j/^ap- As prior function Pr(i;), we initially adopt very 
broad uninformative uniform priors: for the logarithmic slope we 
consider the interval ye [1,3]; for the lens strength cq e [0, 1]; the 
inclination can assume all the values i e [0°, 90°]; finally, the axial 
ratio q 6 [0, 1.5] can account for all the flattenings between a thin 
disk and a very prolate distribution. Fast preliminary exploration 
runs then peimit us to select, for each system, narrower priors that 
remain, however, generously non-restrictive (i.e. they contain the 
bulk of the posterior and a wide margin of safety all around it). 

The uncertainties on each individual parameter rji, are obtained 
by marginalizing the posterior PDF over all other parameters. We 
refer to the parameter value that corresponds to the maximum of 
the A;-th one-dimensional marginalized distribution as qk.mux - In this 
context, it is important to note that, in general, //map do not 

coincide and can actually differ significantly. In other words, the 
MAP parameters individually are not guaranteed to be probable: 
although the MAP solution corresponds by definition to the highest 
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value of the joint posterior PDF, it does not necessarily occupy a 
large volume in the multivariate parameter space and could easily 
lie outside of the bulk of the posterior PDF if the distribution is 
strongly non- symmetric. 

The computationally expensive task of exploring the poste- 
rior distribution Pr(j7 1 d) is accomplished by m aking use of the 
very efficient and robu st MultiNest algorithm jperoz & HobsonI 
l2008l :l Feroz et al. I l2009l) , which impleme nts the 'nested samp ling' 
Monte Cario technique Jskilling 2004 ; Sivia & Skilling|2006f) . and 
can provide reliable posterior inferences even in presence of multi- 
modal and degenerate multivariate distributions. For the analysis of 
the SLACS lenses, we launch MultiNest with 200 live points (the 
live or active points are the initial samples, drawn from the prior 
distribution, from which the posterior exploration is started) to ob- 
tain the individual marginalized posterior PDFs of the non-linear 
parameters ag, y and q. In the following, unless stated otherwise, 
we quote the 99 per cent confidence intervals (calculated such that 
the probability of being below the interval is the same as being 
above it) obtained from these distributions as our uncertainty on 
the parameters. 



3.4 Results 

The CAULDRON code has been applied to the analysis of all six- 
teen lens galaxies for which VIMOS integral-field spectroscopic 
data are available. These objects are all early-type galaxies, with 
the exception o f J 125 1 -0208, which is moiphologically classified 
by lAuger et"al] ( |2009|) as a late-type galaxy. We include the lat- 
ter in the E/SO sample studied in this paper, since it is a bulge- 
dominated system which was spectroscopically selected from the 
SDSS using the same criteria — optimized to detect bright early- 
type lens ga laxies — adopted for the other objects (see lSolton et al.l 
l2005l l2006h . For six of the systems, the results of the combined 
lensing and dynami c s study were reported in previous publications 
dCzoske et Zll2008l : iBamabe et alj|2009al) . However, updated and 
extended kinematic maps are available for J2321 and therefore we 
have re-analyzed that galaxy. 

Table[T]lists, for each galaxy, the recovered non-linear param- 
eters for the MAP model, i.e. the inclination /, the lens strength ao, 
the logarithmic slope y' and the axial ratio q of the total density 
distribution. The reconstructed observables corresponding to this 
model are presented and compared to the data sets in Appendix |A] 
(available in the online version of the journal), where also the re- 
covered unlensed background source is shown; the MAP model 
reconstruction of the weighted DF is presented in Sect. |5] where 
we study the dynamical structure of the lens galaxies. None of the 
galaxies requires external shear to be added to the model. 

The marginalized posterior PDFs for the non-linear model pa- 
rameters — which quantify the statistical uncertainties on those pa- 
rameters (see Sect. 13. 3t — are presented in Appendix lB](available 
in the online version of the journal) for all sixteen galaxies in the 
sampl^. In Table [T] we also indicate, as a compact description of 
the eiTors, the parameter values which correspond to the maximum 



* The error analysis for the galaxies studied in iBamabe et al.l )2009ah was 
performed using a less robust implementation of the nested sampling tech- 
nique. For consistency, the posterior distributions for those systems have 
been re-analyzed here using the more robust MultiNest algorithm. The re- 
vised confidence intervals are found to be wider, as a consequence of the 
improved parameter space exploration. 
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Figure 1. The logarithmic slope of the total density profile plotted against 
redshift for the sixteen eai'ly-type lens galaxies in the ensemble. The error 
bars are calculated by considering the region of the marginalized posterior 
PDF for y' (see Figure IBU which contains 99 per cent of the probability. 
The blue line indicates the slope y' = 2, corresponding to the isothermal 
profile, while the dashed black line denotes the average slope (y') = 2.07 for 
the sample. The grey band region visualizes the intrinsic scatter of about 7 
per cent in y' (see text for details). 



of the marginalized posterior distributions (subscript 'max') and the 
limits of the 99 per cent confidence intervals. 

All the observables are, in general, accurately reproduced 
within the noise level. In particular, the lensed image residual maps 
do not exhibit the prominent structured patterns that can be seen 
when the method is applied to complex simulated systems, as in 
the [Sarnabe et al. (2009b) tests, thus indicating that the underly- 
ing density distribution of real ellipticals is fairly smooth and well- 
behaved, and can be satisfactorily modelled by the axisymmetric 
power-law profile of Eq. ([T}. 

The reconstructed background sources generally appear as 
simple and smooth systems dominated by a single component, but 
a few of them display a more complex morphology or fainter sec- 
ondary components. To understand this, we recall that the SLACS 
candidates are spectroscopically selected from the SDSS database 
by identifying those systems which have composite spectra consist- 
ing of both an absoiption-dominated continuum (due to the fore- 
ground galaxy) and multiple nebular emission lines at a higher 
redshift (due t o the background o bject), which directly trace star 
formation (see lBolton et aljr2006h . Therefore, the background ob- 
jects tend to be actively star-forming galaxies which, in the typical 
redshift range where the sources are located, i.e. 0.4 < z < 0.8, 
can sometimes be late-type systems charac terized by clumpy, ir- 
regular an d flocculent appearances (cf. e.g. IZamoiski et al.ll2007l ; 
lElmegreen et al. 2009). 

The complicated and patchy moiphology of certain sources 
(as in the case of SDSSJ2303 or, more conspicuously , of the 
SLACS lens SDSS J0728 analyzed in lSamabe et ai]|20I(]|) as well 
as the presence of multiple sources almost adjacent to each other 
(e.g. SDSSJ1251, SDSSJ2321) can thus be interpreted as dis- 
tinct star-forming regions within the same galaxy. Analogously, the 
small clumps observed in SDSSJ1250 can likely be identified as 
bright blue star-forming regions that are observed within an inter- 
acting system with a tidal tale. 
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Table 1. Recovered parameters for the sixteen sample SLACS lens galaxies. 



'/MAP 

Galaxy name z\cm Jsrc 



CIO q I y ao 



SDSS J0037-0942 





1955 


0.6322 


1 


968 





434 


0.693 


65 


5 


1 


nf,Q+0.014 





400+0.018 

-^-^-0.009 


0.694+«;«JJ 


65 


4+2 
^-0 


SDSSJ0216-0813 





3317 


0.5235 


1 


973 





344 


0.816 


70 





1 


079+0.035 
' -0.077 





044+0.0 18 
-0.012 


0.818!«;0« 


69 


51 


SDSS J09 12+0029 





1642 


0.3240 


1 


877 





412 


0.672 


87 


8 


1 


ssn+0.035 
^'^"-o.oie 





4,0+0.005 
^'^-0.004 


670+'' °2l 

-0.026 


87 


6- 


SDSS J0935-0003 





3475 


0.4670 


2 


285 





225 


0.982 


34 


8 


2 


900+0.072 
'^'*"-0.056 





990+0.013 
^^°-0.019 


0.982+003| 


34 




SDSSJ0959+0410 





1260 


0.5349 


1 


873 





323 


0.930 


80 


4 


1 


o£Q+0.013 





09,^+0.009 
•'^-'-0.008 


0.920!«;gi 


80 




SDSSJ1204+0358 





1644 


0.6307 


2 


226 





399 


1.000 


65 


2 


2 


990+0.057 
"-0.074 





009+0.020 
^**"^-0.011 


0.986+«;?g« 


65 


8+« 


SDSSJ1250+0523 





2318 


0.7953 


2 


131 





339 


0.771 


29 


1 


2 


, ^4+0.025 
" -0.074 





336+°-''2' 

-0.011 


n 779 +0.054 
"■"'^-0.133 


28 




SDSSJ1251-0208 





2243 


0.7843 


1 


925 





244 


0.796 


85 


6 


1 


□24+0.009 
'^^-0.030 





945:+0.005 
-^-0.006 


n on-^+oois 

-0.046 


85 


4+] 


SDSSJ1330-0148 





0808 


0.7115 


2 


276 





155 


0.414 


74 


3 


2 


9sn+0.039 





1 1:4+0.008 
'-^ -0.012 


0.414+«;0^« 


75 


3- 


SDSSJ1443+0304 





1338 


0.4187 


2 


435 





149 


0.698 


71 


5 


2 


411+0.069 
-0.076 





15Q+0.020 
-0.021 


0.682+«;0« 


71 


8-T 


SDSSJ1451-0239 





1254 


0.5203 


2 


053 





332 


0.999 


40 


8 


2 


090+0.030 
02s 





007+0.003 
•'•"-0.007 


n 997+0.012 
"■'"-0.055 


40 


8- 


SDSSJ1627-0053 





2076 


0.5241 


2 


122 





369 


0.851 


56 


4 


2 


190+0.029 
'^-^-0.016 





O6Q+0.004 
^"'-0.011 


n = 0+0.014 


55 


8+1 


SDSSJ2238-0754 





1371 


0.7126 


2 


088 





362 


0.781 


79 


8 


2 


nQO+0.022 
"^'^-0.053 





oco+O.Oll 
-^•^ -0.005 


767+'' 042 
"■'"'-0.067 


79 


1+' 

-2 


SDSS J2300+0022 





2285 


0.4635 


1 


921 





349 


0.642 


59 


4 


1 


090+0.033 
' -^-0.043 





044+0.007 
•^ -0.004 


622+'' 034 
"■"^^-0.012 


58 


7+4 
-1 


SDSSJ2303+1422 





1553 


0.5170 


2 


102 





436 


0.642 


89 


1 


2 


nQo+0.030 

"y»-0.004 





407+0.003 
^•"-0.007 


644+0 004 
"■"^^-0.033 


85 


9+3 

'-0 


SDSS J232 1-0939 





0819 


0.5324 


2 


058 





467 


0.744 


67 


3 


2 


Q59+0.013 

'-0.066 





4/;o+0.009 
4&**-0.003 


744+0.029 

"■' -0.011 


67 


5+2 

•^-2 



Note. The non-linear r] parameters are: the logarithmic slope 7'; the dimensionless lens strength ao; the axial ratio q and the inclination 1 (in degrees) with 
respect to the line-of-sight. Columns 4 to 7 list the MAP parameters, i.e. the parameters that maximize the joint posterior distribution ('best model' parameters). 
Columns 8 to 11 encapsulate a desciiption of the one-dimensional marginalized posterior PDFs (shown as histograms in Appendix iBlin the online version of 
the journal): rj^^i^^ is the maximum of that distribution and the indicated en^ors represent the lower and upper limits of the 99 per cent confidence interval. 




Figure 2. Joint posterior probability for y[ and CTyi (see Eq. (7)). The cross 
marks the position of the joint maximum. The three contours con^espond 
to the posterior ratios PlPmux = e^^^~^^, with A^^ = 1,4,9. The contours 
are only for indication: formally, they can be inteipreted as 1 , 2 and 3 cr 
contours only for a (multivariate) Gaussian distribution, in which case 
represents the usual chi-square. 



3.4.1 The slope of the density profile 

By examining the recovered values for the logarithmic slope of the 
total density distribution, we find no evidence of evolution of 7' 
with redshift within the probed range (z = 0.08 to 0.35), confirm ing 
the results of iKoopmans et all ( |2009| ) and I Auger et al.l j2010ah on 
the full SLACS sample (see Figure [T). Moreover, from the same 
Figure, we also see that, while almost all systems have a nearly (i.e. 
within 10 per cent) isothermal profile, very few of them exhibit a 
slope that is actually consistent with y' = 2 within the uncertainties. 

We now proceed to study in a more quantitative way the dis- 
tribution of density slopes for the ensemble at hand. Under the as- 
sumption that the slope y[ (with a symmetrized 1 cr error 5y[) of 
each lens system is drawn from a parent Gaussian distribution of 
centre y'^ and standard deviation ay>, the joint posterior probability 
P of these two parameters for the sample is given by 



'P=-Px{y[,ay\[y'^)<xp{y[,ay) 







exp 





(7) 



where p(y'^,(Ty>) indicates the prior, for which we adopt a uni- 
form distributiorlfl This probability function, visualized in Figure|2] 
yields the average logarithmic density slope 



^ Since the results are essentially data driven, the specific choice of the 
prior distribution is largely irrelevant. If, for instance, we adopt a scale- 
invariant Jeft'reys prior p oc I /cry (which indicates the absence of ci priori 
infonnation on the scale of iry ) the recovered MAP value for o-y differs 
less than 4 per cent from the value obtained above with p = constant. 
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Figure 3. Axial ratio q of the total density distribution plotted against the 
intrinsic axial ratio </» of the luminous distribution. For each galaxy, q, is 
calculated by deprojecting the observed isophotal axis ratio 5».2D: using the 
recovered MAP value of the inclination. The dashed line corresponds to 
q = q,. Red and blue symbols denote, respectively, slow and fast rotators. 



</> = 2.074!»:«« (68% CL), 

which is slightly super-isothermal. This result is con s istent with 
the value of {y') d etermined by the iKoopmans eFai] j2009l) and 
lAuger et"al] ( l2010ai) combined lensing and dynamics analysis of 
the full SLACS sample (which does not include two-dimensional 
kinematic information), using spherical Jeans equations and assum- 
ing an isotropic stellar velocity dispersion tensor 

The intrinsic scatter around the average slope, also derived 
from the posterior distribution of Eq. (|7), is 



: 0.144! 



(68% CL), 



which corresponds to 6.9!^^ per cent of (y'), and is displayed 
in Figure [T] as a grey band. This value is only marginally 
lower than the intrinsic sca tter of about 10 per cent obtained by 
iKoopmans etai] jlOOSl) and I Auger et"!!] j2010j) . and the two (Tyi 
are consistent within the errors. 



3.4.2 The flattening of the density profile 

All the analyzed galaxies — with the exception of SDSS JI330 — 
are found to be in the range between moderately flattened and per- 
fectly spherical, with a total density axial ratio 0.6 < q ^ 1.0. 
The ensemble-average axial ratio, obtained from the joint posterior 
probability as described in Sect. |3.4.T| is 



{q} = 0.770! 



(68% CL), 



with a quite substantial intrinsic scatter of about 19 per cent around 
this value: 

= 0.148!o;ot5 (68% CL). 

The three roundest galaxies in the sample, and the only 
ones consistent with being spherical, are all slow rotators (see 
lEmsellem et al.ll2007l and Sect. |5.2| for our definition of slow and 
fast rotating systems based on the specific angular momentum). 



The system SDSS J1330, which is peculiar for its remarkably flat- 
tened axial ratio q =^ 0.4, also bears the distinction of being the 
least massive object and one of the only five fast-rotating galaxies 
in the sample. Fast rotators, however, are not necessarily very flat- 
tened: SDSS J0959, for instance, is found to be almost spherical, 
with q = 0.9. 

In Figure |3] ^ is compared to the intrinsic (i.e. three- 
dimensional) axial ratio of the luminous distributions, q,, calcu- 
lated as 



l-(l-ql,„)/smH, 



(8) 



where q, 2D is the observed isophotal axial ratio and / indicates the 
recovered MAP value for the inclination. For most systems, the 
flattening of the luminous distribution is consistent or very similar 
(i.e. within ~ 10 per cent) to the flattening of the total distribu- 
tion; we note that this does not necessarily imply that mass fol- 
lows light as the total and luminous radial density profiles could 
be very different. Only three galaxies exhibit a significant discrep- 
ancy, all of them in the sense of having a total axial ratio much 
rounder than the luminous one, i.e. q/q, > 1.4. Of these systems, 
two (SDSS J0959 and SDSS J1251) are fast rotators, whereas the 
third one, SDSS J0935, is a very massive slow rotator (and the most 
massive galaxy in the whole sample, cf. Sect.|4ll. Out of the remain- 
ing fast rotating galaxies, SDSSJ1443 shows q/q, =^ 1.2, while 
SDSS J 1330 and SDSS J2238 are slightly rounder in the total distri- 
bution but still consistent with q = q,.In conclusion, from the small 
sample at hand, we find that fast rotators tend to be more flattened 
in the stellar mass distribution than in the total density distribution 
(i.e., q is at least as round as q,, and generally rounder). 



4 MASS STRUCTURE 
4.1 Total mass 

For each galaxy, the (spherically averaged) total mass profile cor- 
responding to the best reconstructed model is calculated from the 
density distribution of Eq. ([TJ, and shown as the solid black curve 
in Figure|4] The vertical lines indicate, for reference, the location of 
the (two-dimensional) effective radius R^, the Einstein radius i?Einst 
and the outermost radius iJ^n probed by the kinematic data set. For 
most systems the kinematic maps extend up to and beyond the half- 
light radius, providing strong constraints within this region, and in 
all cases (with the exception of SDSS J0912) Ri^i„ exceeds at least 
RJ2. 

While /?Einst and /fkin provide useful yardsticks to compare the 
data coverage in difi^erent systems, it should be noted that the ob- 
servational constraints to the mass models a re not strictly limite d 
to these radial values since, as discussed in ICzoske et al.l l l2008h . 
the data sets also include the contribution of more distant galaxy 
regions that are located along the line-of-sight. 

The analyzed sample spans almost two orders of magnitude in 
mass: the total mass Mto,,e enclosed within the three-dimensional 
radius r = R^, hereafter denoted as re for brevity, ranges from a few 
lO'^M© for the two smallest systems SDSS J1330 and SDSS J1443 
to the almost lO'^M© of SDSS J0216 and SDSSJ0935, with most 
of the galaxies falling in the range 1 - 3 x 1O"M0 (see Table|2]l. 

There is some evidence (non-zero slope with about 2-sigma 
significance) for a slight anti-coiTelation between the total mass 
M,o, e and y' (see Figure |5] upper panel, and Table |3), i.e. the to- 
tal density slope tends to be steeper for the least massive systems. 
The most notable outlier is SDSS J0935, which, despite being the 
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most massive galaxy in the sample, has a y' comparable to that 
of the two systems with M,ot,e < 10"Mo. We also find a cor- 
relation (with similar significance) between y' and the quantity 
log(M,o,i./r^), which is a proxy for the average total mass density 
inside the three-dimensional radius re (F igure [5] lower panel) . This 
trend is in agreement with the findings of I Auger et al .1 ( 120 10a ). who 
observe a relatively tight correlation between y' and the central sur- 
face mass density of the full SLACS sample. 

We note that the covariance between y' and M,(,( e is generally 
very small. In fact, for a given Einstein mass and radius, the de- 
pendence on the slope of the total mass Miat(r) enclosed inside the 
three-dimensional radius r is least in the region where r/REimt - 2. 
For most of the galaxies in our sample, the ratio i^c/^Einst is not 
far from that value, so that the 1-sigma uncertainty on y' usually 
translates into a mass change of less than 1 per cent (and even in 
the case of J0935, which has ife/^Einst - 5, it only reaches about 4 
per cent, i.e. less than 0.02 in dex). 



4.2 Luminous mass and dark matter content 

The radial profile for the luminous component is obtained from 
the best reconstructed stellar DF. However, since the adopted mod- 
elling approach treats the stellar distribution as a tracer embedded 
in the total potential, the normalization of the luminous profile re- 
mains unconstrained, and further assumptions (§ 14.2.1b or addi- 
tional information (§ 14.2.21 ) are therefore necessary to set its scale 
and determine the contribution of the dark matter component at dif- 
ferent radii. 



4.2.1 Dark matter fraction: constraints from the maximum bulge 
approach 

When additional information on the normalization is unavailable, 
the most straightforward approach consists in maximally rescaling 
the luminous profile under the physical constraint that it does not 
exceed the total mass distribution at any point, and the assumption 
that the stellar mass-to-light ratio is fairly position-independent. 
This procedure is known as the minimum halo or maximum bulge 
hypothesis, and is the analogue for early-type galaxies of the clas- 
sical maximum disk hypothesis l lvan Albada & Sancisilll98^ rou- 
tinely adopted in t he study of late-type gal axies. It is frequently 
employed (see e.g. iBamabe et al ] |2009al and lWeiimans et al.ll2009l 
for recent applications) since it provides a consistent and robust 
method to determine a conservative lower limit for the dark matter 
fraction in the galaxy, within the considered region. 

The (spherically averaged) luminous mass profiles obtained 
under the maximum bulge hypothesis are presented as dashed black 
lines in Figure |4] Let M™'^ denote the corresponding stellar mass 
enclosed inside the spherical radius re. We see that the lower limit 
for the (three-dimensional) fraction of dark over total mass within 
that radius, defined as /dm = 1 - M™e''/Mtot,e, varies very signif- 
icantly from system to system (cf. Table |2}: three of the galax- 
ies, SDSS J1443, SDSS 11204 and SDSS J2303, show a dark matter 
lower limit /dm < 5 per cent within the probed region, whereas — 
on the opposite end — several galaxies require at least one fifth of 
the total mass enclosed within r^ to be dark, with a peak of 46 per 
cent for SDSS J1251. The average and median of the dark matter 
fraction lower limits for the sample, obtained under the maximum 
bulge approach, are, respectively, /dm = 16 and 12 per cent. 

The stellar mass-to-light ratios (in the 6-band) corresponding 
to the maximal luminous profiles are in the range 3 < M./Lq s < 



i 2 




log[M,„^,/10'"M3,J 



i 2 




log [ (M 



Figure S. Logarithmic slope of the total density profile y' against the total 
mass enclosed within the three-dimensional radius (e (upper panel) and the 
average total mass density inside re (lower panel). The solid black line is a 
hnear fit to the relation (including scatter), and the gray band indicates the 
intrinsic scatter 



10, fully consistent with the values determined for local early- 
type galaxies (see e.g. Kronawitter et a l. 2000 Gerha rd et alj|200ll 
iTmjillo. Burkert. & Bellll2004 iMamon & LokasI I2005h . However, 
when taking into account passive stellar evolution, galaxies at 
higher redshift are expected to have smaller M,/LQ g than their 
local count erparts (an effec t that becomes significant already at 
z > 0.1, see lTreu et al.l2002h . 



4.2.2 Dark matter fraction: constraints from stellar mass 
estimates 

While calculating a lower limit for the dark matter content is use- 
ful, galaxies need not have maximal bulges. We therefore set the 
scale for the luminous mass profiles by adopting, for each system, 
the stellar mass in ferred from stell ar po pulation synthes is models, 
assuming either a lChabried (l2003h or a ISalpeterl (Il955l) IMF. We 
use stellar mass values taken from the stellar population analy- 
sis performed by lAuger et al.1 ( |2009|) on a data set constituted by 
deep, high-resolution, multi-band HST observations of the SLACS 
lenses. 

The spherically averaged luminous mass profiles obtained in 
this way are presented in Figure |4] Figure |6]shows the correspond- 
ing dark matter fractions within r^, following the definition for /dm 
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Figure 4. Compilation of spherically averaged mass distributions for the sample galaxies, sorted by increasing Mtoi e . In each panel, the solid black line shows 
the total mass profile obtained from the best reconstructed model, and the dashed black line indicates the maximal luminous mass profile calculated under 
the maximum bulge (i.e. minimum h alo) hypothesis (se e Sect. I4.2."n . The solid red line shows the luminous mass profile rescaled using the value for stellar 
mass determined from SPS analysis jAuger et al,ll2009l) . adopting a Salpeter IMF, The red-shaded region represents the corresponding uncertainty (i.e., the 
standard deviation of the marginalized posterior for the stellar mass). The blue curves and regions have the same meaning, but are obtained with a Chabrier 
IMF. The vertical lines indicate the location of the effective radius (dotted fine), the Einstein radius (dashed line) and the outermost boundary of the kinematic 
data (dash-dotted line). 



given above in § 14.2.11 and replacing Mf J" with the stellar masses 
M'^^'^ and A/f e'' for Chabrier and Salpeter IMFs, respectively. 

We find that, independently of the specific choice of the IMF, 
there is a significant correlation (greater than 3-sigma) between the 
dark matter fraction within re and the total mass (the parameters 
of the fit are given in Table [3ll. Such a correlation was first ob- 
served, in the SLACS sample, bv lAuger et al. I( l2010al) . Their result 
is confirmed and strengthened by the present study, where we con- 



duct a more rigorous and detailed joint analysis on a sub-sample 
of SLACS galaxies, taking advantage of the extended kinematic 
information (rather than an estimate of the stellar velocity disper- 
sion (TsDss within a single 3 arcsec diameter aperture), and self- 
consistently including axial symmetry and stellar anisotropy in the 
model, and we reach analogous conclusions. 

The normalization based on the stellar masses inferred from 
Salpeter IMF produces a luminous profile that, when compared to 
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Table 2. Recovered masses, dark matter fractions and dynamical quantities for the sixteen sample SLACS lens galaxies. 



Galaxy name ll|!f ^'I^Q.B /dm tim foM "I"' h ^ 



SDSS J0037-0942 


11.50 


4.85 


SDSS J0216-0813 


11.97 


9.80 


SDSS J0912+0029 


11.91 


9.00 


CT^CC THOQ^ 


1 1.97 


7.17 


SDSS J0959+0410 


10.99 


7.63 


SDSSJ1204+0358 


11.14 


7.55 


SDSSJ1250+0523 


11.29 


3.55 


SDSSJ1251-0208 


11.27 


3.56 


SDSSJ1330-0148 


10.29 


5.55 


SDSSJ1443+0304 


10.62 


4.04 


SDSSJ1451-0239 


11.11 


4.04 


SDSSJ1627-0053 


11.37 


4.85 


SDSS J2238-0754 


11.16 


5.91 


SDSS J2300+0022 


11.44 


7.22 


SDSSJ2303+1422 


11.51 


7.48 


SDSSJ2321-0939 


11.27 


4.81 



0.23 


0.66+«o^ 


40+008 

-0 09 


0.17 


0.83;«»3 


0.70;0-M 


0.30 


0.77!«M 


0.59!0-06 


0.12 


0.81!«;03 


67+0-05 


0.30 


-0.05 


49+0-0' 


0.02 


57+" °'^ 

"■-^ -0.08 


23+0-'0 


0.05 


Q 35+010 


-0 14+0- " 


0.46 


67+'' 05 

' -0.07 


43+0.10 

"■^^-0.11 


0.05 


51+0.06 

"■■^ -0.07 


15+0-" 

■ -0.12 


0.01 


32+0 09 


-0 20+0"' 

"■^"-0.18 


0.12 


51+0.07 

-0.09 


19+0-'0 

■ '-0.12 


0.21 


50+0 09 


12+0-"' 


0.08 


62+0 05 


n :i^+0.09 
"■^^-0.10 


0.23 


69+0 05 

^'■"^-0.05 


45+0-08 

"■^-"-0.10 


0.04 


69+0 04 


46+0-0' 

"■^"-0.08 


0.12 


51+0.08 


13+0- '5 

"■'^-0.18 



0.08 


0.25 


0.05 


0.16 


-0.37 


0.04 


0.12 


0.02 


0.08 


-0.17 


0.12 


0.23 


0.13 


0.07 


-0.15 


0. 10 


0. 17 


0.03 


0.24 


—0.64 


0.49 


0.65 


0.42 


-0.16 


0.27 


0.03 


0.07 


0.03 


0.09 


-0.19 


0.15 


0.16 


0.09 


0.22 


-0.56 


0.82 


0.81 


0.66 


-0.36 


0.53 


0.52 


0.49 


0.48 


-0.05 


0.09 


0.36 


0.60 


0.37 


0.13 


-0.30 


0.10 


0.31 


0.05 


0.15 


-0.35 


0.07 


0.18 


0.06 


0.16 


-0.38 


0.43 


0.66 


0.42 


-0.18 


0.30 


0.03 


0.15 


0.02 


0.08 


-0.17 


0.05 


0.19 


0.04 


0.14 


-0.31 


0.06 


0.08 


0.06 


0.15 


-0.36 



Note. For each galaxy we list: the logarithm of the total mass Mmt.c enclosed within the three-dimensional radius ; the stellar mass-to-light ratio (in the 
i?-band) corresponding to the maximal luminous profile; the fraction of dark over total mass within re for the maximum bulge assumption (/^jjj ) and for 
the Chabrier and Salpeter IMFs {f^^ and , respectively); the inclination-corrected «/(T ratio; the dimensionless angular momentum j-; the Emsellem 
kinematic parameter Ar; the global anisotropy parameters S and y. 



Chabrier, is always closer to the one determined under the maxi- 
mum bulge assumption, and in several cases is fully consistent with 
it0 For two systems, SDSS J1250 and in particular SDSS J1443, 
this luminous profile (including the error band) locally exceeds the 
maximum bulge curve, indicating that it unphysically overshoots 
the total mass in places, thus producing the two negative /dm points 
in Figure |6] Taken at face value, this would lead to the conclu- 
sion that the Salpeter relation constitutes an inadequate IMF for 
these two galaxies (and therefore must be rejected for all systems if 
we assume the IMF to be universal), whereas the Chabrier IMF — 
which gives lower stellar masses by almost a factor of 2 — never 
fails to produce physical results. However, it must be noted that in 
both cases the lower limit of the error band is only slightly higher 
(< 5 per cent) than the maximal lum inous curve, and th e adopted 
uncertainties represent l-cr errors (see I Auger et al .120091 for the de- 
tails of how these uncertainties are calculated within a Bayesian 
framework). Moreover, several systematic effects (such as the cir- 
cularization of the profile and the discreteness of the TIC building 
blocks used to construct the model galaxy) might interfere at the 
few per cent level with the minute details of the reconstructed lu- 
minous profile. Therefore, we conservatively conclude that all the 
analyzed systems are consistent with having stellar masses deter- 
mined from a Salpeter-like IMfQ 



'' The SLACS lens galaxy SDSSJ0728 (analyzed in lBamabe erai]|201(]|) . 
which has kinematic maps obtained from Keck spectroscopy, represents, 
however, an exception to this trend, in the sense that for that galaxy the 
luminous mass profile obtained with a Chabrier IMF almost coincides with 
the maximum bulge determined one. We note that this systems lies close to 
the low end of the probed range in velocity dispersion. 
' The situation, i.e. a recovered luminous profile slightly overshooting the 



With a Salpeter IMF, the recovered dark matter fraction for 
systems having total masses log(Miot,c/MQ) < 11.5 spans from 
to ~ 50 per cent, whereas the three most massive galaxies are al- 
ready dominated by dark matter (/dm is approximately 60 to 70 per 
cent) inside 1 r^. The average and median dark matter fractions are, 
respectively, = 31 and 37 per cent. The importance of the dark 
component becomes even more extreme if a Chabrier IMF is con- 
sidered (average /p,^* = 61 per cent, median f^^^ = 64 per cent), 
with the luminous component contributing less than 50 per cent to 
Mjot.c in all systems but two. 

It has been suggested that the IMF is not univer- 
sal but becomes sy stematicall y 'hea vi er' for mo re mas- 
sive galaxies (see iTYeu et alj l2009l lAuger et al. ] l2010bl 



Ivan Dokkum & Conrov 2010l 



Thomas et 



lAuge: 



201 1 ; but see 



also Napolitano. Romanowskv, & TortoraT I2OIOI for a different 
conclusion). The data presented here are consistent with this fact 
and cannot resolve the degeneracy between varying IMF and 
varying dark matter fraction, as we discuss in Section |6] 



5 DYNAMICAL STRUCTURE 

A major strength of the method described in Sect.|3]is that the caul- 
dron algorithm also provides, for each analyzed galaxy, the associ- 
ated (weighted) two-integral stellar DF, which constitutes the most 
general and complete description of the dynamical structure of the 
system within the context of the adopted model. The DF can equiv- 
alently be represented as a bi-dimensional map of the TIC weights 



physical limit of the maximum bulge, and the conclusions are analogous for 
the Keck system SDSS J0728 mentioned above. 
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Figure 7. Best model reconstruction of the weighted two-integral DFs of the sample galaxies. The value of each pixel represents the relative contribution of 
the coiTesponding TIC to the luminous mass of the galaxy. The galaxies are sorted by increasing values of the specific angular' momentum parameter j-. 



in the integral space (E, L-): the value of each pixel simply indi- 
cates the relative stellar-mass contribution due to the TIC building 
block of corresponding energy and angular momentum. 

Figure |7] shows the maps of the TIC weights that are ob- 
tained for the best-reconstructed model (i.e., the MAP model) of 
each lens galaxy. In keeping with the analysis of SDSSJ0728 
( iBamabe et allllOIOl) . the most recent study of a lens galaxy per- 
formed with CAULDRON, We make use of a library composed of 360 
TICS: for each one of the A'^ = 20 elements logarithmically sam- 
pled in the circular radius (corresponding to a grid in energy, see 



ISamabe & KooDmansl2007h we consider N^, = 9 elements linearly 
sampled in the normalized angular momentum ^ L./L,^^^^ ^ I, 
and the mirrored negative L-IL-^,^ values. Although a grid of 
{Ne = 10) X (Nl. = 5) has been shown to be generally suffi- 
cient to provide a satisfactory reconstruction of the observables 
iBarnabe et al. 2009a), a finer grid does a better job in reducing 
the undesired discreteness effect of the TIC superposition (e.g. on 
the surface brightness, kinematic and local {v^)l& maps). 

The DF encapsulates the description of the stellar component 
of the galaxy in an very compact way, which can not be straightfor- 
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Figure 6. Dark matter fractions versus the total mass enclosed within 
the three-dimensional radius r^. The upwards black triangles indicate the 
dark matter fraction lower limits calculated from the maximum bulge ap- 
proach. The full circles show the dark matter fractions obtained when lu- 
minous masses are determined via SPS analysis, assuming Salpeter (red) or 
Chabrier (blue) IMFs. The dashed lines are linear fits to the relation, calcu- 
lated including intrinsic scatter (which is indicated by the coloured bands). 



Table 3. Summary of linear relations derived for the SLACS lenses. 



X 


Y 


Slope 


Intercept 


Scatter 


log(Mtoi,e) 

log(Mto,e/r3) 

log(M,o,.e) 
log(M,o,.e) 


y' 
Y 

/Chab 
•'DM 
^S:ilp 
•'DM 


-0.16 ±0.09 
0.16 ±0.09 
0.22 ± 0.04 
0.39 ± 0.08 


2.29 ±0.13 
1.94 ±0.08 
0.35 ± 0.06 
-0.13 ±0.12 


-0.01 
14+0 06 

''■'^-o.oi 
05+0 04 

10+0 06 

"■'"-0.04 



Note. The mass Mmt.c is in units of 10'"Mo, while the density Mtot,c/'e is 
in units of 10^ Mq kpc"^. 



wardly decoded at a glance. Therefore, we now proceed to distill 
much of that information into more directly interpreted quantities 
of astrophysical interest that characterize the dynamical status of 
the system on a global and local level, i.e. the ratio between or- 
dered and random motions, the specific angular momentum and the 
global anisotropy of the stellar velocity dispersion tensor. 

5.1 Global and local u/cr 

Figure [8] shows the location of the sixteen analyzed lens galaxies 
on the model version of the classical (u/cr, e) diagram. The quantity 
vja is commonly used as a global indicator of the importance of 
ordered motions (i.e. rotation) with respect to random motions in a 
stellar system. For each galaxy we calculate the ratio via, corrected 
to an edge-on view, from the two-dimensional kinematic maps of 
the be st reconstructed model, following the prescriptions of lBinnevI 
( l2005h for extended kinematic data sets. The ratio is plotted against 



the inclination-corrected ellipticity of the light distribution 6, = 1 - 
qt, where q, is the axial ratio introduced in Sect. 13.4.2] 

We find that the dynamics of most systems is clearly domi- 
nated by random motions, with via ^ 0.15, while rotation plays 
a significant role in the remaining five galaxies. The members of 
the latter group have several common characteristics, beyond being 
above a via ratio of 0.35: (i) they are all fast rotators, according to 
both the definition of ISmsellem et alj | |2007|) and their specific an- 
gular momentum content (as discussed below in Sect. l5.2t ; (ii) they 
are quite flattened in the light (having 0.38 < q, < 0.74); (iii) they 
belong to the lower half of the sample in terms of Mto,_e, includ- 
ing, moreover, the three least massive galaxies of the whole sample 
(namely, SDSS J1330, SDSS J1443 and SDSS0959). 

As a comparison to the SLACS systems. Figure |8] also dis- 
plays the location on the diagram of a sub-sample of 24 SAURON 
early-type galaxies, whose observables ha ve also been corrected for 
inclination, analyzed by ICappellari et al] ( l2007b . The two samples 
generally occupy similar positions in the (via, e) space, although 
the SLACS galaxies lack the most extreme cases (such as the very 
flattened slow rotator NGC4550 and the fast rotator NGC3156 
with a via ~ 1) and show a sharper transition between dispersion- 
and rotation-dominated systems. This, however, is likely a con- 
sequence of the limited number of fast rotating galaxies present 
among the SLACS sample, mainly comprised of massive systems 
(o"sDss It 200 km s"') which tend to be slow rotators. 

From the stellar DF it is also possible to derive a description of 
the intrinsic kinematics of the galaxy by using the internal velocity 
moments to define — at each point {R, z) in the meridional plane — 
the quantity {v^)la, where {v^) is the local mean azimuthal velocity 
(i.e. around the symmetry axis) and a^ = {a\ + a^ + cr?)/3 is the 
mean velocity dispersion. This ratio can be interpreted as a local 
and intrinsic analogue of the global indicator via, since it charac- 
terizes the importance of rotation with respect to random motions at 
each position inside the galaxy. Figure|9]shows the (v^)la maps for 
all the sample galaxies, extending up to r^ll and sorted by increas- 
ing specific angular momentum /; (see Sect. 15. 2t . For visualization 
purposes, and in particular to make the comparison among fast ro- 
tators easier, the maps are presented as having the same sense of 
rotation (where applicable). 

By examining the {v^)la maps, the sample galaxies clearly 
fall into two broad groups on the basis of their dynamical status. 
The first 1 1 systems are overall dominated by random motions; in 
a few systems, such as SDSSJ0935 and SDSSJ1451, islands of 
moderate rotation {\{v^)\la a 0.5) are present, but remain confined 
to small, spatially limited regions. What sets apart the remaining 
five galaxies, on the other hand, is the presence of strong rotation 
on a large scale pattern. This structural difference reflects the clas- 
sification of the two groups as, respectively, slow and fast rotators. 
Guided by the intuition provided by these maps, in the next Sec- 
tion we provide a more quantitative criterion, based on the angular 
momentum content, to specify to which one of the two groups the 
modelled galaxies belong. 



5.2 The angular momentum of slow and fast rotators 

We recall the definition of the dimensio nless specific angular mo- 
mentum of the stellar component j- (see lBamabe et al.ll2009ah : 

r p,R\(v^)\A'x 

r I (9) 

I p,R^{v^f +a^d^x 
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Figure 9. Maps of the local (v^)/<t ratio between the mean rotation velocity around the i-axis and the mean velocity dispersion, plotted up to Rc/2 in the 
positive quadrant of the meridional plane. As in Fig.|7] the galaxies are sorted by increasing j-. 



where p, = f fi^v is the density of the luminous component ob- 
tained from the MAP model stellar DF /. For each galaxy in the 
sample, we carry out this integral up to re/2, tabulating the results 
in Table |2] We find that the galaxies can be neatly divided in two 
groups depending on the value of the quantity j-: on one side 11 
low angular momentum systems, with j- < 0.3, and on the other 
side five objects with j, 0.5 or greater By comparing with the 
local {Vf)l& maps, we see that these five systems are also the only 
ones that clearly exhibit large scale rotation. Therefore, we choose 
to classify the galaxies on the basis of their specific angular mo- 
mentum content by labelling as slow and fast rotators the systems 
that are, respectively, below and above j- = 0.4. 



In Figure [TO] we plot, for each system, j- versus the quan- 
tities via (inclination-corrected. Sect. I5.lt and Ar, computed on 
the kinematic data set of the best reconstructed model. The ob- 
servational parameter Ar was first introduced by femsel lem et al.l 
( i2007i) as a consistent way to quantify the specific stellar angu- 
lar momentum (in projection) of an early-type galaxy for which 
two-dimensional kinematic maps are available, and its value is em- 
ployed as the criterion to discriminate between slow (Ar < 0. 1) and 
fast (Ar > 0.1) rotating systems. The model quantity /. is, by con- 
struction, the three-dimensional and intrinsic analogue of Ar. Fig- 
ure [TO] shows that classifying the rotators on the basis of either of 
these two parameters yields fully equivalent results for all galaxies 
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Figure 8. Model (u/cr, e) diagram for the sixteen lens galaxies in our sam- 
ple (open circles), divided in fast and slow rotators (blue and red points, 
respectively); e, is the intrinsic elhpticity of luminous distribution, and u/cr 
is calculated from the best model, connected to an edge-on inclination. As 
a comparison, the plot also shows, as full circles, the corresponding quan- 
tities (als o con'ected for inclinatio n) for the 24 nearby SAURON ellipticals 
studied in ICappellari erai]i2007l) . The solid hne r epresents the location of 
edge-on isotropic rotators (assuming o- = 0. 15, see lBinnevll2005l) . 



by 7 = 26/(6 - 1). The latter scenario is the one that applies to 
the galaxy models considered in this analysis, since, for any ax- 
isymmetric coUisionless system supported by a two-integral DF, 
{vj,) = (or) everywhere, the meridional velocities (vr) and {v-} are 
null, and therefore Hhr = n~. 

The values of the global parameters 6 and y for the SLACS 
galaxies, calculated within are reported in Table |2] The slow 
rotators are all mildly anisotro pic, quite similarly to w hat is found 
for their nearby counterparts (ICappellari et al.lcOO?!) : the values 
of 6 fall between 0.05 and 0.25, and 6 < 0.15 for most systems. 
The fast rotating galaxies, with the exception of SDSS J1443, ex- 
hibit instead negative 6 values. This is a consequence of the adopted 
dynamical model: due to the isotropy in the meridional plane, it fol- 
lows from Eq. dlOt that every system for which cr^ < cr^ over most 
of the luminosity density- weighted volume (so that H^^ < Tim) is 
characterized by 6 < 0, i.e. is anisotropic in the sense of having 
a smaller pressure along the equatorial plane than perpendicular 
to it. Moreover — within the approach of modelling based on TIC 
superposition — very fast rotating systems tend to have small az- 
imuthal velocity dispersions cr^, since their dynamical description 
is dominated by co-rotating TIC building blocks of similar angular 
momentum. This can be more easily illustrated by considering the 
limiting case of a hypothetical system described by a single TIC: 
for this simple object one has (d^) = (t)^)^, and thus cr^ = 0, ev- 
erywhere in the meridional plane within the zero-velocity curve, 
so that 6tic = -1 and yxic = 1. It is not surprising, therefore, 
that the three fastest rotating galaxies in the sample (SDSS J0959, 
SDSS J2238 and SDSS J1251) are all found to have cleariy negative 
values of 6. 



but one, SDSS J0912. We have already singled out this system (in 
Sect.llll as the only one for which the kinematic maps do not reach 
half of the effective radius. The misclassification of SDSSJ0912 
can be ascribed to the spatially limited data set, since slow rotators 
commonly show decreasing Ar profiles for R/R^ < 0.5, and thus 
the paramet er can exceed the thr eshold value of 0.1 in the inner 
regions (see lEmsellem et al .120071 and in particular their Figure 2). 

It is worth noting that, for slow rotators, the spread in Af. is 
much smaller than the spread in indicating that the latter con- 
stitutes a better discriminator of the intrinsic properties of these 
systems. 



5.3 Global velocity dispersion tensor 

The global stellar velocity dispersion tensor constitutes a useful 
and concise characterization of the overall dynamical structure of 
a stellar system. It is customary to descri be its shape by means 
of the three global anisot ropy parameters dCappellari et al. I [20071 ; 
iBinnev & Tremainell2008l) 



1 - 



and 6=1- 



2n„ 



n™ -I- n,„ 



(10) 



where we denote the total unordered kinetic energy along the coor- 
dinate direction k as 



(11) 



and crfix) = {v\)-{Vk)^ is the local stellar velocity dispersion in the 
k direction. For a fully isotropic object the values of all three pa- 
rameters become zero. If the galaxy is only isotropic in the merid- 
ional plane, then /? = and the remaining parameters are related 



6 DISCUSSION 

The study presented in this paper represents the most detailed anal- 
ysis conducted so far of the total mass density profile of a relatively 
large sample of E/SO galaxies beyond the local Universe, making 
consistent use of axially symmetric density distributions for both 
the lensing and dynamical modelling, and taking advantage of the 
constraints provided by two-dimensional kinematic information. 



6.1 Total density profile 

The results of this work corroborate the mounting evidence (see 
Sect.lTJ that the inner regions (on the scale of the effective radius) of 
massive early-type galaxies are characterized by an approximately 
isothermal total density distribution (i.e., ptot °^ > with y' » 2), 
a property often referred to as the 'bulge-halo conspiracy' . This 
name owes to the fact that there appears to be no obvious physical 
explanation of why the luminous and dark matter components of 
each system, despite not obeying a power-law distribution, should 
often 'conspire' to produce a nearly l/r^ total profile. 

By examining individual lens galaxies in great detail, how- 
ever, we can conclude that the purported conspiracy only holds in 
first approximation: even though we provide confirmation that all 
systems, quite remarkably, can be effectively modelled by a simple 
and smooth single power-law profile for ptot, the average logarith- 
mic slope is found to be (y') = 2.074;^j5[j^j, i.e. slightly (but signifi- 
cantly) super-isothermal. Moreover, and perhaps more importantly, 
different galaxies do exhibit noticeably different slopes, with an in- 
trinsic scatter in (y') of ab out 7 per cent. These re sults are in agree- 
ment w ith the findings of iKoopmans et all i200Sl ) and lAuger et al.l 
l l2010al) for the full sample of SLACS lenses, obtained by making 
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Figure 10. The intrinsic dynamical parameter j- versus the inclination- 
corrected v/cT ratio (top panel) and the Emsellem kinematic parameter Afi 
(bottom panel) for the sixteen SLACS lens galaxies in the sample. The cir- 
cles and diamonds indicat e, respectively, the slow and fast rotating systems 
classified according to the lEmsellem et ai]i2007l) definition of the threshold 
value that divides the two types, i.e. = 0. 1 (horizontal dashed line). The 
red and blue symbols denote the slow and fast rotators when the galaxies are 
classified on the basis of the dimensionless angular momentum (threshold 
value = 0.4, vertical dotted line). The two classifications produce fully 
equivalent results for all systems barring SDSSJ0912 (indicated with the 
red diamond), which is the galaxy with the most spatially limited kinematic 
maps (see text). 



use of spherical Jeans models for the dynamics, constrained by a 
single kinematic measurement, i.e. the stellar velocity dispersion 
(TsDss inside the 3 arcsec diameter SDSS fiber. This indicates that 
the combined lensing and dynamics approach is a generally robust 
method to determine the total density distribution of a sample of 
lens galaxies, even when a very simple (and not fully consistent) 
dynamical model is employed. The values of the logarithmic slope 
obtained for individual systems with the two methods — that we 
denote here as y^i and y'j^.^^^ — can differ by up to about 10-15 per 
cent, as illustrated in Figure [TT] This discrepancy is typically of 
the order of the 1-cr errors in yj^.^,,, but, in most cases, much larger 
than the few per cent uncertainties on y' obtained with this study, 
stressing the importance of having access to two-dimensional kine- 
matic data sets and more sophisticated dynamical models when a 
precise measurement of a lens galaxy density profile is called for. It 
is interesting to note that another joint lensing and dynamics study 
conducted bv lRuff^ et al. (20 11) on a different sample of 11 massive 
early-type galaxies from the Strong Lenses in the Legacy Survey 
(SL2S) at a higher redshift (median zjens = 0.5) also finds an av- 
erage slope and an intrinsic scatter consistent with the results pre- 
sented here. 




Figure 11. Comparison between the total density logarithmic slope of 
the sample lens galaxies obtained from the analysis conducted in this pa- 
per, making use of two-integral dynamical models a nd VLT VIMOS inte - 
gral field spectroscopy, and the slope yjj.^^^ derived bv lAuger et al.l l llOlOall . 
using spherical Jeans models and a single kinematic measurement. Red 
and blue points denote, respectively, slow and fast rotators. The dashed 
line traces the one-to-one relationship. For both y^j and y';^.^^^^, the uncer- 
tainties represent the 68 per cent confidence interval. The slope of galaxy 
SDSS 1251 is not listed, as it was not included in the I Auger et al.l i2010al) 
sample of early-type systems. 



It has been well established by N-body cosmological simu- 
lations that dissipationless processes produce inner density distri- 



butions with logarithmi c slopes y' x I (e.g. iNavarro et al 



or even shallower (e.g. iGraham et al.l 12006 : iNavarro et al 



1996b 



20101) . 



Moreover, subsequent dissipationless merging (i.e. , dry merging) 
preserves the steepness of the inner density profile (lDehnenll2005l : 
iKazantzi dis et al . ,200d). It follows, therefore, that the interplay 
between dark matter and the baryonic component, with its com- 
plex dissipative processes, plays a key role and must be taken 
into account in order to explain the emergence of the isothermal 
profile in the inner regions of galaxie s. Numerical work (see e.g. 
iGnedin et aLll2004 lAbadi et al.ll2010l and references therein) on 
the response of dark matter to the condensation of baryons does 
indeed show that the assembly of a central galaxy makes the halo 
significantly more concentrated, producing a steeper total density 
slope. The high-resolution simulations of moderately massive sys- 
tems (with halo masses compa rable to that of the Milky Way) car- 
ried out bv lTissera et al. l (l2O10h . which include a detailed treatment 
of baryonic processes, reveal a nearly isothermal profile over a rel- 
atively large radial range extending up to the baryonic radius (de- 
fined as the radius containing 83 per cent of the stellar and gaseous 
component of the galaxy, i.e. the region within which the baryons 
play an important role). Aside from this general global result, how- 
ever, the details of the density distribution are largely determined 
by the specific assembly history of each system. This suggests 
that the observed scatter in the y' values of the SLACS lenses is 
a consequence of the idiosyncratic formation histories of the in- 
dividual galaxies, and can potentially provide information on the 
processes involved in their build-up. Once the isothermal struc- 
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ture is in place, dissipationless mergers — both major and minor — 
reserve it not only in tlie inne r regions, but throughout the systems 
Nipoti. Treu. & Boltonl2009h . We stress, however, that, even in the 
case of a nearly perfectly 1 /r^ seed distribution, this process pro- 
duces a non-negligible scatter of about 10 per cent in the logarith- 
mic slope (i.e. comparable to the observed level), which overlaps 
with the scatter that, potentially, might have been present in the 
parent (pre- dry mergers) density distribution, further complicating 
the task of reconstructing the galaxy formation history. 



6.2 Mass and dynamical structure 

Another defining property of early-type galaxies that we have 
aimed to quantify in this study is the fractional amount of dark 
matter that is present in the systems' central regions. If the assump- 
tion of a constant stellar mass-to-light ratio holds, at least within 
the inner or so, then most of the analyzed systems do require 
a non-negligible fraction of dark matter (/dm greater than 10 per 
cent) even when the contribution of the luminous component is 
maximized ad hoc, without any concerns for the plausibility of the 
MJLqs ratios involved in the rescaling. In other words, what is 
shown by this 'maximum bulge' approach is that the mass-foUows- 
light hypothesis is necessarily violated in the majority of cases. 

Adopting a physically motivated rescaling of the luminous 
profile — based on the stellar masses inferred from SPS models — 
leads to interesting results. 

Fi rst of all, in agreement with previous work o n the SLACS 
lenses jKoopmans et alj llOOd : iBolton et all l2008hl ; I Auger et all 
l2010al : lBamabe et al.ll201oir we find a significant amount of dark 
matter within the galaxy inner regions: in the case of a Salpeter 
IMF, the average /dm is 31 per cent, with a wide variation 
within individual objects, ranging from galaxies consistent with 
having no dark matter, to systems where the baryonic compo- 
nent accounts for less than half of the total mass already at r^. 
These values are broadly consistent with the results of earlier dy- 
namical studies of nearby early-type galaxie s employing a vari' 



ety of different technioues and assumptions jGerhard et al 
Cappellari et al.ll2006l ; lThomas et al.ll2007l : IWeiimans et~ 



2001 



Tortoraetal.ll2009l) . which find dark matter fractions inside 



200S 



half-light radius varying from to about 50 per cent of the to- 
tal mass, with a typical average value close to 30 per cent. If we 
prescribe, in stead, an IMF lighter th an Salpeter — as advocated for 
example bv ICappellari et al]l2006l based on their analysis of the 
SAURON ellipticals, which have lower velocity dispersions than 
the SLACS lenses considered here — then the dark component be- 
comes the dominant one in the central regions for almost all sys- 
tems (average /dm is 61 per cent with a Chabrier IMF), in strik- 
ing contrast with the traditional picture of early-type galaxies. It is 
worth noting that similarly abundant dark matter contents within 
the galaxy inner regions are commonly f ound by studi es of E/SO 
galaxies that adopt a Chabrier-like IMF: iGrillol ( |2010|) . e.g., cal- 
culates a median value of 64 per cent (in projection) within 1 R,. 
for a large sample o f 1.7 x 10' massive SDSS early-type galaxies. 
iTortora et all | |2009|) . applying spherical Jeans models to a homo- 
gene ous sample of 33 5 local E/SO systems, derive a /dm ~ 60 per 
cent. IWeiimifflj | [2009l) performs a sophisticated dynamical analysis 
of the nearby galaxy NGC 2549 using high-quality photometry and 
extended kinematic maps and finds, likewise, a dark matter fraction 
of 65 per cent within the half-light radius by adopting a iKroupal 
( I2OOII) IMF 

The second relevant result is the non-negligible positive cor- 
relation between the derived dark matter fraction and the total mass 



enc losed within rg. This relatio n is in agreement with the results of 
the IPadmanabhan et al] | |2004|) study of a sa mple of about 30,000 
SD SS ellipticals, and with the findings of iTortora et al.l 
and lOraves & FabeJ ( 1201 Ol) at lower redshift. Conversely, 

— by determining the stellar masses from SDSS multicolor 
photometry and the total mass from a simple dynamical relation, 
under the assumption that the galaxies can be described as spheri- 
cal and isothermal systems — comes to a different conclusion, i.e. 
that the (projected) dark matter fraction within the half-light radius 
remains almost constant for values of the total mass between a few 
I0'° and ~ 10'- Mq. 

If the IMF is truly universal, then the trend for /dm implies a 
genuine increase of the dark matter contribution with M,oi.e, and for 
the most massive ellipticals the stars represent a minority compo- 
nent, even in the central regions (cf. Sect. 14.2.2) . However, an al- 
ternative explanation is possible if the IMF varies with galaxy mass 
from a Chabrier/K roupa-l ike functional form to a Salpeter-like one 
jTreu et al]l201(ll : see also lAuger et al.ll2010bl). or even to a s teeper 
low-mass end slope, as argued bv Ivan Dokkum & Conrovl l l2O10h 
for the most massive systems (i.e., cr > 250 km s"', which would 
include many of the objects studied here). We note, however, that 
the results from the SLACS survey only imply a mass-to-light ra- 
tio that is higher than expected for a Chabrier or Kroupa IMF. The 
excess mass can be provided by lo w mass stars (as in the case of a 
Salpeter IMF, as advocated by, e.g.. lvan Dokkum & Conrovll2010l) . 
or by high mass stars remnants, in th e form of neutron star and 
black holes jTreu et al.ll2010l : lAuger et al . 2010b). 

Finally, in terms of dynamical structure, the SLACS lenses 
present no surprises, and are found to be very si milar to their lo- 
cal counterparts as observed , e.g., by SAURON dEmsellem et al.l 
I2OO7I : ICappellari et al.ll2007l) . More than two-thirds of the galax- 
ies are slow-rotating objects, as a consequence of the sample be- 
ing skewed towards the high velocity dispersion end. The five fast- 
rotators largely coincide with the low mass tail of the sample, and 
all of them are high angular momentum systems, with Ar well 
above 0.1, and very clearly identified as such when their {v^}/d- 
meridional plane maps are examined. The dark matter fraction of 
the fast-rotators is consistent with the general trend, and they do 
no t appear as a different p opulation in terms of /dm. as suggested 
bv lCappellari et al.l ( 12009) . 

There is tentative evidence that fast-rotating galaxies are sys- 
tematically more flattened in the luminous distribution than in the 
total density profile, despite the fact that these systems — as dis- 
cussed above — are among the least massive in the sample and thus 
likely to be the ones where the baryonic component is most impor- 
tant. This would then imply that these objects lie in dark haloes 
whose central regions are significantly rounder than the light we 
observe. Such a scenario appears to be in general agreement with 
the results of nu merical simulations (e.g. iDebattista et al.l I2OO8I : 
iAbadi et al.ll2010h , which show that the assembly of massive high 
angular momentum galaxies changes the dark halo into nearly ax- 
ially symmetric, slightly oblate systems over a wide radial range. 
The question of whether there is an actual structural difference be- 
tween the dark haloes of slow and fast rotators, however, remains 
pending. 



7 SUMMARY 

In this work we have conducted a combined, self-consistent grav- 
itational lensing and stellar dynamics analysis of the full sample 
of SLACS lens systems with available VLT VIMOS integral-field 
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spectroscopic data, employing axially symmetric dynamical mod- 
els supported by two-integral DFs. This is integrated with the con- 
straints from stellar mass estimates derived from SPS analysis of 
multiband HST imaging. The application of these modelling tools 
on this remarkable data set has enabled us to investigate in de- 
tail the properties and three-dimensional structure of the consid- 
ered sample of sixteen early-type galaxies in the redshift range 
z = 0.08 - 0.35. 

This study represents an improvement over previous joint 
lensing and dynamics analyses (which adopt dynamical models 
based on Jeans equations) and constitutes, in a sense, an extension 
beyond the local Universe of SAURON-type studies of the inner re- 
gions of E/SO galaxies. Although the detailed kinematic maps avail- 
able for nearby systems are clearly superior to the corresponding 
data sets currently obtainable for objects at z > 0.1, the additional 
information derived from gravitational lensing provides an invalu- 
able asset for the study of distant galaxies, which allows us to put 
robust constraints on their mass distribution and amount of dark 
matter within the probed regions. 

The main conclusions from this analysis are summarized as 
follows. 

(i) The total mass density distribution of all the sample galaxies 
is well described, within the inner regions probed by the data sets 
(of order one half-light radius), by an axially symmetric power-law 
model. The average logarithmic slope is slightly super-isothermal, 
with (y') = 2.074+Jj j5;|J; there is an intrinsic scatter cry = 0.144+55-555J 
(corresponding to less than 10 per cent) around this value. 

(ii) We find that the lens galaxies have a fairly round total mass 
distribution; the average axial ratio is (q) = 0.77 ± 0.04 and — 
while there are significant differences between individual systems, 
with an intrinsic scatter of nearly 20 per cent — all objects, barring 
one, are rounder than q = 0.6. Most galaxies are about as flattened 
in the total distribution as they are in the luminous one. 

(iii) The lower limit for the fraction of dark over total mass 
within the three-dimensional radius re (calculated with the maxi- 
mum bulge approach under the assumption of constant stellar mass- 
to-light ratio inside the considered region) varies significantly be- 
tween individual systems from nearly zero to almost 50 per cent, 
with an average value /dm = 16 per cent (median value: /dm = 12 
percent). The M,/Lq^ ratios corresponding to these maximal light 
profiles range from 3 to almost 10. 

(iv) When the normalization of the luminous profile is set using 
stellar mass determinations obtained from SPS models, under the 
assumption of a universal IMF of the Salpeter type, we derive an 
average dark matter fraction f^'^ = 31 per cent inside r^ (median 
value: f^'^ = 37 per cent). None of the analyzed galaxies has a 
luminous mass distribution strongly at odds with the assumption of 
a Salpeter IMF (the luminous profile does exceed unphysically the 
circularized total mass profile for two systems, but in both cases the 
discrepancy is just a few per cent above the l-cr uncertainty). 

If one adopts a Chabrier IMF the recovered dark matter frac- 
tion becomes much higher, i.e. /qm* = 61 per cent (median value: 
fmf' ~ ^ P^'^ cent), which would imply that the stars represent a 
minority mass component even in the inner regions of early-type 
galaxies. 

(v) If the stellar masses are determined in this way, i.e. assuming 
a universal IMF, then the derived /dm shows a clear coiTelation 
(3-sigma) with M,at,c, in the sense that the dark matter contribution 
increases with the total mass of the galaxy, becoming — even in the 
case of Salpeter IMF — the dominant component within r^ for the 
most massive ellipticals (i.e. the systems with log[Mto,_i:/M0] > 



11.5). If the IMF is not universal, than a progressive steepening of 
its profile with the galaxy total mass could explain, at least partially, 
the observed trend. 

(vi) The SLACS lenses can be divided into two dynamically 
distinct groups based on the value of their specific angular mo- 
mentum parameter j- (which we can introduce based on our three- 
dimensional models): about two- thirds of the galaxies are domi- 
nated by random motions and show little large-scale ordered ro- 
tation (0 < /; < 0.3), whereas the remaining five, all of which 
are among the least massive systems in the sample, are high an- 
gular momentum objects (/ ^ 0.5 and above). These two groups 
correspond to the standard classes of slow and fast rotators, as re- 
vealed by several other properties (e.g. mass range, flattening, the 
Emsellem parameter Ar), and in particular by the characteristics of 
the respective {v^}/d- maps, which illustrate the local distribution of 
ordered-to-random motions ratios in the meridional plane. For slow 
rotators, the spread in / is much larger than the spread in Ar, sug- 
gesting that the specific angular momentum parameter is a better 
discriminator of the intrinsic properties of these systems. 

(vii) Our analysis shows that the SLACS lens systems are over- 
all analogous to their local Universe counterparts in terms of den- 
sity profile and global structural and dynamical properties. This in 
turn indicates that massive early-type galaxies have experienced at 
most limited structural evolution, at least as far as their inner re- 
gions are concerned, within the last four billion years. 
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APPENDIX A: BEST MODEL RECONSTRUCTIONS FOR 
THE SAMPLE GALAXIES 

In this Appendix we present thie data sets (i.e., galaxy-subtracted 
lensed image, surface brighitness maps and kinematics maps), tlie 
corresponding best model reconstructed observables, and the resid- 
uals for eleven systems in the considered sample. The correspond- 
ing panels for the remaining five galaxies (namely J0037, J0216, 
J0912, J0959 and J 1627) were pres ented in a previous publication 
in this series jBamabe et al. I l2009a|) and therefore will not be du- 
pli cated here. The res ults for galaxy J2321 were first described 
in jCzoske et al.ll2008h ; however, since an improved analysis of 
the J2321 kinematic data set has become available, we have re- 
analyzed that system and we have included it in the selection pre- 
sented here. 

A thorough description of the data sets will be provided in a 
forthcoming publication (Czoske et al. 201 1, in preparation). 
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Figure Al. Best model lens image reconstruction for the system 
SDSSJ0935. From the top left-hand to bottom right-hand panel: recon- 
stracted source model; HST/ ACS data showing the lens image after subtrac- 
tion of the lens galaxy; lens image reconstruction; residuals. In the panels, 
North is up and East is to the right. 
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Figure A2. Best dynamical model for the galaxy SDSS J0935. First row: 
observed surface brightness distribution, projected line-of-sight velocity 
and hne-of-sight velocity dispersion. Second row: conesponding recon- 
structed quantities for the best model. Third row: residuals. In the panels. 
North is up and East is to the right. 
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Figure A3. Best model lens image reconstruction for the 
SDSS J1204. Panels meaning as in Fig.lAT] 
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Figure A4. Best dynamical model for the galaxy SDSS J 1204. Panels mean- 
ing as in Fig. lA2l 
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Figure A5. Best model lens image reconstruction for the galaxy 
SDSS J1250. Panels meaning as in Fig.lAT] 
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Figure A6. Best dynamical model for the galaxy SDSS J 1250. Panels mean- 
ing as in Fig. lA2l 
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Figure A7. Best model lens image reconstruction for the galaxy Figure A9. Best model lens image reconstruction for the galaxy 

SDSS J1251. Panels meaning as in Fig.|AT] SDSS J1330. Panels meaning as in Fig.lAT] 
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Figure A8. Best dynamical model for the galaxy SDSSJ1251. Panels mean- 
ing as in Fig. lA2l 



Figure AlO. Best dynamical model for the galaxy SDSSJ1330. Panels 
meaning as in Fig. lA2l 
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Figure A12. Best dynamical model for the galaxy SDSSJ1443. Panels 
meaning as in Fig. lA2l 



Figure A14. Best dynamical model for the galaxy SDSSJ1451. Panels 
meaning as in Fig. lA2l 
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Figure A15. Best model lens image reconstraction for the 
SDSS J2238. Panels meaning as in Fig. lAll 
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Figure A17. Best model lens image reconstiuction for the galaxy 
SDSS J2300. Panels meaning as in Fig.lAT] 
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Figure A19. Best model lens image reconstraction for the galaxy 
SDSS J2303. Panels meaning as in Fig. lAll 
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Figure A21. Best model lens image reconstraction for the galaxy 
SDSSJ2321. Panels meaning as in Fig.lAT] 
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Figure A20. Best dynamical model for the galaxy SDSSJ2303. Panels 
meaning as in Fig. lA2l 



Figure A22. Best dynamical model for the galaxy SDSSJ2321. Panels 
meaning as in Fig. lA2l 
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APPENDIX B: UNCERTAINTIES 

In this Appendix we present, for all the analyzed systems, the one- 
dimensional marginalized posterior PDFs for the considered model 
parameters (i, 7', aa, q), from which the corresponding uncertain- 
ties (quoted in the Table [TJ are derived. 
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Figure Bl. Marginalized posterior probability distribution functions of the power-law model parameters / (inclination), y' (logarithmic slope), ag (lens 
strength) and q (axial ratio) for each of the analyzed systems, obtained from the nested sampling evidence exploration (see text. Sect. (33). From left to right, 
top to bottom, the galaxy are presented in the same order as in Tables[T]and[2] 



